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Abstract 

Relativistic covariance requires that in the constituent quark model for 
mesons the positive energy states as weU as the negative energy states are 
included. Using relativistic quasi-potential equations the contribution of the 
negative energy states is studied for the light and charmonium mesons. It 
is found that these states change the meson mass spectrum significantly but 
leave its global structure untouched. 
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I. INTRODUCTION 



The application of the constituent quark model to mesons and baryons has been very 
successful, in spite of various nonrelativistic approximations. How much is neglected in 
these approximations? Generalization of the ordinary Schrodinger equation to a relativistic 
covariant form results in the well-known Bethe-Salpeter equation (BSE), which differs from 
the nonrelativistic equation in two respects. Firstly, the dependence of the bound state wave 
function on the relative three-momentum p becomes a dependence on the four-momentum 
{Pq,p). Secondly, covariance requires that for fermions the full Dirac structure is taken into 
account, so for quarks not only the positive energy states are to be considered, but also 
their negative energy states must be included. In this paper we investigate the effect on 
the meson mass spectrum of the inclusion of these negative energy states in the relativistic 
constituent quark model. 

Quasi-potential (QP) approximations to the BSE exist which eliminate the po- 
dependence of the quark-propagators, and moreover circumvent the problem of introducing 
a confining potential in the full Minkowski space. The covariant Dirac structure is kept in 
these. In previous work [0,0, to which we shall refer as I and II, we used two QP models 
to calculate the full meson mass spectrum. One of the biggest differences found with the 
Schrodinger approach was the much stronger confinement needed and the sensitivity of the 
mass spectrum to small vector contributions to the confining potential. To get some more 
insight about the various aspects we have performed calculations in the same framework but 
with leaving out various negative energy states. This enables us to trace back and under- 
stand the origin of the phenomena found in I and II and to establish the importance of the 
negative energy states. For simplicity the calculations are restricted to the equal mass case. 

In the next section the model is briefly summarized and some of the QP equations 
properties which give upper bounds on the parameters of the potential are discussed. Section 
III deals with the truncation of the relativistic QP equations in the coordinate represention 
with regards to the negative energy spinor states. It also presents the comparison of the 
calculations with corresponding nonrelativistic predictions. In particular the uu and cc 
spectra are considered in various approximations to various QP models. Differences in the 
spectra are discussed. Finally, the paper closes with some concluding remarks. 



We briefly describe the relativistic constituent quark model as was studied in I and 
II. The model consists of two point-like fermions bound together by a phenomenological 
potential V to form a meson. The meson wave function ip satisfies 



The dependence of the propagator S on the relative four-momentum {pq,p) has been sim- 
plified to a dependence on the three-momentum. In I we studied two QP approximations 
which eliminate the dependence of the equation on pq, the Blankenbecler-Sugar-Logunov- 
Tavkhelidze (BSLT) approximation , and an equal-time (ET) approximation . They 
give for quarks of equal mass in the center of mass system 



II. REVIEW OF MODEL 
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where u = y/p^~+~mF, and E = M/2, M being the total meson mass. The AP'^p^ project 
upon positive and negative energy states, A''^^^ = Ai^{p)A2^{—p) with 
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Let us define the eigenstates of these projection operators by 
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The subscript pw is used to distinguish these plane wave states from the canonical positive 
and negative states defined through 



(1) (2)l \ I \ 

70 70 \Plp2)ca = Plp2\PlP2)ca, 



(2.6) 



which correspond to the various combinations of upper and lower components of the Dirac 
spinors. Clearly 



\PlP2)pw = \PlP2)ca + 0{p/m). 



(2.7) 



The instantaneous interaction V between the quarks is modeled as the sum of a Coulomb- 
like part describing the one-gluon-exchange (OGE) interaction and a linearly rising part for 
the confinement. It takes in coordinate space the form 
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The vector contribution to the interaction is studied in the Feynman gauge as well as in the 
Coulomb gauge: 
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This defines the relativistic quasi-potential model. Note that wc do allow for a fraction e of 
vector- confinement. The running coupling constant behaves as a{x) ~ (8/27)7r/ ln(a;o/a;) for 
small distances x, and grows to some maximum saturation value asat for large separations, 
according to the interpolation given in I. 

The resulting wave equation were studied extensively in coordinate space. It should 
be noted that certain difficulties as found in the one particle Dirac equation should be 
expected in such a relativistic quasi-potential approach. Let us for a moment consider one 
fermion in an external potential. If this particle experiences a potential which fiuctuates 
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more strongly than ~ 2m over a distance shorter than its Compton length xc = 1/m new 
fermion-antifermion pairs can be created. This phenomenon cannot correctly be described 
by the Dirac equation which describes a one-particle theory and thus misses the interactions 
between the newly created pair and the starting particle. Since the Dirac equation does 
allow for antifermion components, solutions in this potential can have an unbound number 
of non-interacting fermions and antifermions thus being unnormalizable and unphysical. 
This break-down of the Dirac equation is well-known as the Klein paradox [0. Similar 
flaws emerge in the QP equations of this work since they also contain negative energy 
components. In view of the complexity of the various two-body equations we do not fully 
analyze under what conditions they break down. Instead, let us note that unbound solutions 
can be expected if the confining strength becomes too strong, kxc ~ 2m or k ^ 2m^. This 
domain can be reached in light meson systems. The condition on k depends on the fraction 
e of vector confinement as indicated by the discussion in I. Similarly, if the OGE potential 
becomes too strong, a/xc <^ 2m or a > 2, irregular solutions may be expected. Mesons with 
high orbital angular momenta are less sensitive to this effect because the centrifugal barrier 
prevents them from entering the short distance region. If a running coupling constant is 
taken instead of a fixed one then the irregular solutions disappear and the singular behavior 
becomes less. Whereas the upper bound on k, disappears as the negative energy components 
are removed, the upper bound on a is shifted upward but still present, reflecting that only 
positive energy states can also tumble in a OGE potential. Detailed discussions on the short 
distance behavior was given in I and II. 



III. RESULTS 



Since the wave equations as described in the previous section are solved in the coordinate 
space using the representation of the canonical states Eq. (|2.6[), a convenient way to switch 
off the coupling to the negative energy plane wave states is to explicitly project out these 
states in the interaction. The projection is done by rewriting the BSLT equation as 
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For = -ibb = 1 we recover the QP equations with all negative energy states included. If i\, is 
put to zero then the single negative energy states are projected to zero, provided E 0, and 

similarly for when we want to drop the coupling to the | )pw states in the calculations. 

One can rewrite the ET equation in a similar way: 



A^ 



E)A++ - uj{A+- + A-+) + (zbbcc' + ^)A— ] ip 
- + Zb(A+- + A-+) + t,i,A—] 7S'^f V^. 



(3.2) 



Although the number of coupled channels do not change, the above Eqs. (|3.1| ) and ( p.2|) 
have a major advantage that they can be solved in the same way as the full BSLT and 
ET equations were solved in I. First we Fourier transform them to configuration space and 
make an angular momentum decomposition. This gives a set of coupled integral-differential 
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equations which are reduced to a set of hnear equations by expanding the wave function on 
a set of sphne functions. The resulting matrix equation for the various sphne coefficients 
can be solved straightforwardly by standard methods. As a check on the accuracy of our 
calculational procedures we used this method to calculate the charmonium mass spectrum in 

the model of Hirano et al. 0, which includes only the |++)p«, and | )pw states. Agreement 

was found within 1 MeV. 



A. Comparison with Schrodinger equation results 



To study the effects of various relativistic contributions we consider the mass spectra for the 
light mesons and the charmonium system. In the limit of large quark masses the BSLT can 
be reduced to the nonrelativistic Schrodinger equation 
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where Vnr is the nonrelativistic reduction of the potential 

Vnr[x) = h + C, 
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and VsD contains the spin-dependent corrections of order 
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Here Vy and Vs denote the vector and scalar contributions to the potential. The potential 
in Eq. ( |3.5| ) is singular at the origin and therefore a regularization is needed. The delta 
function appearing in V^VV is replaced by a gaussian of width /3f ^, and the x~^- and x~^- 
singularities are cut off at xq = and replaced by smoothly fitting gaussians. We take the 
same parameters for the potential as in I, which are summarized in table 1, and take the /3j 
such that they reproduce approximately the same splittings between the S- and P-states as 
the full BSLT calculation. With these parameters we find the spectra as shown in the first 
column of Figs. 1 and 2 of the uu and cc systems. The spectra of the various columns in 
these figures will be referred to as Fig. la, lb, ... etc. 

The first natural step in relativizing the Schrodinger equation is to replace the kinetic 
energy according to 
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in Eq. ( ^.31 ). Figs, lb and 2b show the corresponding spectra. They are shifted downwards 
as compared to the nonrelativistic ones indicating that the interaction has become more 
attractive. Also the shift becomes larger as the level of excitation increases. For charmonium 
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one finds shifts between 0.05 - 0.13 GeV, wliereas for the light mesons these shifts are more 
substantial and range between 0.5 - 1.5 GeV. 

The second relativistic correction which can be included is to replace the approximation 



of the potential Vatj? + Vsd by the complete projection of the potential Eq. ( p.8|) on positive 
energy states V^'^ = pw{++ |^7o^''7o^'' I ++)pw The singular behavior of Vsd is no longer 
present in the complete projection, and hence the cut-off parameters are principly absent 
in the resulting wave equations. Note the factors 7o^''7o V^"*^'"*"*" and in the right hand 
sides of Eqs. ( |3.1| ) and ( p.2| ) which reflect that the energies of the quarks are the fourth 
components of their four-momenta. This replacement gives a rather large effect (0.00 - 0.16 
GeV for cc and 0.0 - 0.8 GeV for uii) as is illustrated by the spectra in Fig. Ic and 2c and 
which is easily understood. For a potential V{x) = Vs{x)I + Vv{x)'yQ'^'y^'^ and low momenta 
of the in- and outgoing states 

V^'-^ ~ ca( + + |^57j'^7f ^ + Vv{x)I\ ++),a = Vs{x) + Vv{x). (3.7) 

But for relativistic momenta considerable contributions to V~^'~^ are to be expected from 
the contributions 

ca{-+ |^57S'^f ^ + Vv{x)I\ -+),a = -Vs{x) + Vv{x). (3.8) 

For scalar confinement it leads to a considerable reduction of the confinement strength. This 
effect can also qualitatively be seen from the spin- independent corrections of order 1/m^ 
to a nonrelativistic scalar potential. It is given by 



Vsiix) = 



(3.9) 



This potential, however, does not give an accurate approximation since for large distances 
mr'^Vs^'^ is not a small enough parameter to expand in. 

To obtain the explicit form of the BSLT propagator we may replace 

{uj~E)^^^{uj~E). (3.10) 
uj + E 

This leads to the BSLT equation restricted to only positive energy states. Figs. Id and 2d 
show how the masses are increased. 

Figure le shows the effect of the introduction of the single negative energy states |H — )pw 
and I — [-)pw For the light mesons we do not show this case since for some mesons it leads 
to unbound systems, similar to the Klein paradox mentioned in the previous section. The 
introduction of these states leads to an increase of the cc spectrum by up to 0.10 GeV. They 
are mostly made out of |H — )ca and | — |-)ca states which have a negative expectation value of 
scalar confinements, Eq. (|3.8|) . Together with the negative propagation of these states this 
leads to a positive interaction which raises the mass levels. 

Finally, if the | )pyj states are included, one arrives at the full BSLT spectra shown 

in Figs. If and 2e. The inclusion of these states changes little (^ 1 MeV for cc), as could 
be expected from their smallness. The lowest lying mesons, however, are considerably in- 
fluenced what is related to the singular behavior of the wave function at small distances. 
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The OGE potential concentrates the wave function to small distances with strong potential, 
thus lowering the mass levels. The singular behavior becomes stronger if more negative 
energy states are included or if the level of excitation is less. Fig. Ig shows the charmonium 

spectrum in the Salpeter approximation, which includes only the |++)p^ and | states. 

The small masses for the low lying excitations are still present, but the increase associated 
with the single negative states is absent. 



B. Equal-time results 

Let us briefly discuss what effects are found if we use the equal time (ET) propagator. 
In Fig. 3 is shown the corresponding charmonium spectrum calculated with the same pa- 
rameters as Fig. 1. Three differences can be noted. Firstly, the ET spectrum is lowered as 
compared to the BSLT one. This is a consequence of the different propagators for the |++)pu; 
states [cf. Eq. ( ^.10[ )]. Secondly, the introduction of the single negative energy states causes 



a shift upwards in the ET model which is approximately twice that seen in BSLT. This is 
due to the fact that these states are roughly twice as important in the ET wave function as 
in the BSLT one in view of a twice as large propagator for these states, S^^ = —l/2u versus 

^BSLT = —l/4:U. Thirdly, the | energy states lower the singlet states much more 

strongly in the ET model than in BSLT. Again, this is a consequence of the importance of 
these states in the ET wave function. They are larger due to the large 5*^^ propagator. As 
already remarked in 0, this is a defect of the ET equation when an unretarded potential is 
used, and this disappears if the proper retardation is inserted in the matrix elements for V 
that connect to double negative energy states. However, it should be stressed that there is 
no unambigious extension of the instantaneous confining potential to retarded times. 



C. Sensitivity to vector confinement 

In 1 we discussed the large raising of the light meson spectrum when a part of vector 
confinement is added to a scalar confining potential. Fig. 4 shows the uu spectrum for 
various dynamical models using scalar confinement without and with a vector contribution. 
When only positive energy states are taken the raising is even more stronger, and can almost 
double some meson masses. 

In Fig. 4 are also shown the same spectrum if the vector contribution is taken in the 
Feynman gauge, Eq. In this gauge the spin-spin interaction is less suppressed so the 

vector contribution has more effect. 



IV. CONCLUDING REMARKS 

In summary, we have calculated for the light and charmonium systems mass spectra 
using the (relativized) Schrodinger equation and the relativistic quasi-potential equations 
of I and II, which differ from the Schrodinger approach in that no approximation is 

made, and that the QP equations contain the full Dirac structure of positive and negative 
energy states. We studied the importance of these differences by solving the QP equations 
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while leaving out negative energy components. We find that the projection of the confining 
potential on positive energy states leads to a considerably lower confinement strength than 
the nonrelativistic potential gives. This is partly compensated for by the introduction of the 
single negative energy states which are more bound than the positive states, and hence in- 
crease the masses. Finally we find that the double negative energy states have httle infiuence 
except for short distances where the singular behavior of the mesons is strengthened. 

The total picture stongly suggests that the differences here studied are important for 
determining the parameters of the qq interaction, especially the confinement strength k. 
Yet the global structure of the spectrum — level ordening, relative sizes of splittings — re- 
mains rather untouched under the relativistic modifications. This confirms the conclusion 
drawn from the success of nonrelativistic quark models that most relativistic effects in qq 
spectroscopy can be mimicked by employing nonrelativistic dynamics together with effective 
parameters. It will be interesting to study whether the relativistic modifications to meson 
wave functions will lead to sizeable changes in cross section for processes such as e.m. ones 
involving these mesons. 

In this paper we have studied relativistic effects within the framework of quasi-potential 
equations. In so doing we have not addressed entirely the role of the relative energy variable 
Pq. Apart from the complexity of a calculation including this say in a Bethe-Salpeter equation 
approach, a more fundamental obstacle is posed by the extension of the definition of the 
confining potential, as we have used here, to a four-momentum dependence. The confining 
potential is only known for the static case. There is at this moment no underlying theory 
which can give a prescription on how to extend it to a covariant form. Let us illustrate this 
considering the commonly used generalization of the potential V{x) = kx which reads in 
momentum space 



92 1 



drf — ql + rf{'^ — is) 



(4.1) 



but its Fourier transform yields 



Vit, x) = - lim \KQir]R) - r]RKAr]R)] = 00, (4.2) 

TT r?-^0 

(with = — t^) which is clearly physically unacceptable. Exploring QCD may lead to 
ways of reconstructing such a confining force, which can be used in such an off mass shell 
approach as we have discussed here. 



ACKNOWLEDGMENTS 

This work was partially financially supported by de Stichting voor Fundamenteel Onder- 
zoek der Materie (FOM), which is sponsored by the Nederlandse Organisatie voor Weten- 
schappelijk Onderzoek (NWO). 



8 



REFERENCES 



[1] P.C. Tiemeijer and J.A. Tjon, Utrecht preprint THU-92/31, |nucl-th 9211003 . 
[2] P.C. Tiemeijer and J.A. Tjon, Phys. Lett. B 277 (1992) 38. 

[3] R. Blankenbecler and R. Sugar, Phys. Rev. 142 (1966) 1051; A.A. Logunov and A.N. 

Tavkhehdze, Nuovo Cimento 29 (1963) 380. 
[4] E.D. Cooper and B.K. Jennings, NucL Phys. A 500 (1989) 553. 

[5] V.B. Mandelzweig and S.J. Wallace, Phys. Lett. B 197 (1987) 469; S.J. Wallace and V.B. 

Mandelzweig, Nucl Phys A 503 (1989) 673. 
[6] S.J. Wallace in Nuclear and Particle Physics on the Light Cone, Proceedings of the 

LAMPF Workshop, Los Alamos 1988. 
[7] see e.g. J.D. Bjorken and S.D. Drell, Relativistic quantum mechanics (McGraw-Hill, New 

York, 1964). 

[8] M. Hirano, K. Iwata, K. Kato, T. Murota and D. Tsuruda, Prog. Theor. Phys. 69 (1983) 
1498. 

[9] W. Lucha, F.F. Schoberl and D. Gromes, Phys. Rep. 200 (1991) 127. 



9 



FIGURES 



5.5 



S=0, J=L 



5.0 - 



4.5 



> 
O 



S=1,J=L+1 

S=l, J=L 

S=1,J=L-1 



S 4.0 - ----- 



3.5 - 



3.0 




a 



f 



g 



FIG. 1. Charmonium spectrum of the radially unexcited states, a: Schrodinger equation 
with p^/2m and Breit interaction, b: Schrodinger equation with -y/p^Tjl^^ ^nd Breit interaction, 

c: Schrodinger equation with ^/p"^ + w? and full projection of potential in Coulomb-gauge, i.e. FT 



with I ++; 



d: BSLT with 



+ +)pw, e: BSLT with 



+ +)pw, |H — )pw and | — h) 



ipwi 



f: BSLT with 



all states, g: BSLT with | H — \-)pw and 



Ipw 
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FIG. 2. Light meson spectrum of the radially unexcited states, a: Schrodinger equation 
with /2m and Breit interaction, b: Schrodinger equation with \Jp^ + m? and Breit interaction, 
c: Schrodinger equation with ^/p^ + m? and full projection of potential in Coulomb-gauge, i.e. ET 
with I ++)pw, d: BSLT with | ++)pw, e: BSLT with all states. 
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FIG. 3. Charmonium spectrum of the radially unexcited states (legend as in Fig. Ill) 



a: Schrodinger equation with /2m and Breit interaction, b: Schrodinger equation with \/p'^ + m? 
and Breit interaction, c: Schrodinger equation with ^/p^ -\- m? and full projection of potential, i.e. 

+ +)pw, |H — )pw and I — l-)piu, e: ET with all states, f: ET with 



ET with I d: ET with 



+ +)pw and 



Ipw 
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FIG. 4. Light meson spectrum of the radially unexcited states (legend as in Fig. |^). a: BSLT 
in Coulomb gauge with | ++)pw and e = 0, b: BSLT in Coulomb gauge with | ++)pw and e = 0.25, 
c: BSLT in Coulomb gauge with all states and e = 0, d: BSLT in Coulomb gauge with all states 
and e = 0.25, e: BSLT in Feynman gauge with | ++)pw and e = 0, f: BSLT in Feynman gauge with 
I ++)pw and e = 0.25. 
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TABLES 



TABLE L Parameters for the qq model 





uu 


cc 




m 


0.250 


1.779 


GeV 


K 


0.33 


0.33 


GeV^ 


c 


-LO 


-1.0 


GeV 


£ 


0.25 


0.25 




asat"" 


0.8 


0.8 




Pi 


0.25 


0.8 


GeV 




0.10 


0.6 


GeV 


^Running as aj in Ref. | 


I with fj, = 1.0 
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